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^ ' Abstract. In the spin- i frustrated chain with nearest-neighbor ferromagnetic exchange Ji and 

C/3 , next-nearest-neighbor antiferromagnetic exchange J2 under magnetic field, magnetic multipolar- 

hquid (quadrupolar, octupolar, and hexadecapolar) phases are widely expanded from the 
saturation down to a low-field regime. Recently, we have clarified characteristic temperature 
and field dependence of the NMR relaxation rate l/Ti in the quadrupolar phase. In this 
paper, we examine those of 1/Ti in the octupolar phase combining field theoretical method 
with numerical data. The relevance of the results to quasi one-dimensional J1-J2 magnets such 
as PbCuS04(OH)2, Rb2Cu2Mo30i2 and Li2ZrCu04 is shortly discussed. 
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1. Introduction 



(N 
> 

^ . 1 . 

, The spin- 2 chain with nearest-neighbor exchange couphng Ji and next-nearest-neighbor couphng 

' J2 is one of the representatives of frustrated spin systems [T]. If J2 is antiferromagnetic (AF) 

i.e., J2 > 0, the system has magnetic frustration regardless of the sign of Ji. Recently, several 
quasi one-dimensional (ID) edge-sharing cuprates [2l[3llll[5l[6l[7l[8l[9] such as LiCu202, 
LiCuV04, Rb2Cu2Mo30i2, PbCuS04(OH)2 and Li2ZrCu04, have been investigated as quantum 
multiferroic materials or frustrated magnets, and many of them are believed to be approximately 
described by a J1-J2 spin chain with ferromagnetic (FM) Ji and AF J2 [IHl [11] at least in the 
high temperature regime above their ordering temperatures. 

Recent theoretical studies have elucidated effects of external magnetic field H on the J1-J2 
' chain [12^ I13j. It was found that the field induces three multipolar Tomonaga-Luttinger (TL) 

liquid phases, i.e., quadrupolar, octupolar, and hexadecapolar TL-liquid phases, in the wide 
region of {Ji/ J2, H). In these three phases, multipolar operators nn=i ^J+n '^ith p = 2, 3, and 
4 are respectively shown to exhibit a quasi long-range order. In addition, in these phases, the 
longitudinal spin correlator also decays in an algebraic form, while the transverse spin correlator 
obeys an exponential-decay law. 

It is generally difficult to experimentally distinguish magnetic multipolar phases in spin 
systems. There is no established method of measuring multipolar correlations, and static 
quantities such as specific heat and susceptibility exhibit no large difference between the 
multipolar and standard TL liquids. Recently, we have shown [2] for these multipolar TL liquids 
that (i) in contrast to usual TL-liquid phase, nuclear-magnetic-resonance (NMR) relaxation rate 



1/Ti of the multipolar liquid phases decreases with lowering temperature near the saturation, 
and (ii) gapless points and peak positions of dynamical structure factors are clearly different 
from those of usual AF spin chains. Furthermore, in Ref. [15] we have numerically evaluated 
1/Ti of the quadrupolar phase as a function of temperature and magnetic field, combining 
field theoretical approach with numerical data. The obtained result demonstrated that field 
dependence of l/Ti is also characteristic: 1/Ti decreases with increasing the field except for 
the vicinity of saturation. These features of dynamical quantities would be helpful to detect 
multipolar TL liquids in experiments for the above quasi-lD magnets. In this paper, using the 
same method as Ref. [TJ], we determine the field and temperature dependence of 1/Ti in the 
octupolar phase of spin-| J1-J2 chain. The result is expected to be relevant for quasi-lD J1-J2 
magnets with large FM |Ji|. For instance, PbCuS04(OH)2 might be a possible candidate for 
the realization of octupolar TL liquid |16] . 

2. Properties of multipolar phases 

Here we shortly review the low-energy properties of multipolar TL-liquid phases of spin-^ J1-J2 
chain [12[ 113]. The Hamiltonian of J1-J2 chain under magnetic field H is defined as 

H= ^ ^J„,S,.S,+„,-^FS|, (1) 

n=l,2 j j 

where we focus on the case with Ji < and J2 > 0. For the saturated state in high- field 
regime, two-, three- and four-magnon bound states with momentum vr, respectively, become 
the lowest-energy excitation in the regions < IJ1I/J2 < 2.699 (0.371 < J2/IJ1I < 00), 
2.720 < IJ1I/J2 < 3.514 (0.285 < J2/IJ1I < 0.368), and 3.514 < IJ1I/J2 < 3.764 (0.266 < 
J2/\Ji\ < 0.285). With decreasing field H below the saturation field, these bound magnons 
are condensed, and the system can be regarded as a hard-core Bose gas of bound magnons. 
We call these gases of two-, three-, and four-magnon bound states, respectively, the second, 
third, and fourth order multipolar (quadrupolar, octupolar, and hexadecapolar) TL liquids [T7] . 
These multipolar phases expand from the saturation down to the low-field regime and turn into 
the low-field vector-chiral phase with {{Sj x Sj^iY) 7^ 0. In the p-th order multipolar phase, 
correlation functions of both the longitudinal spin operator Sj and the p-th multipolar operators 
]^^^-^ "^j+n (i-^-' creation/annihilation operator of bound magnons) can be calculated from the 
effective TL-liquid theory for the hard-core Bose gas of multi-magnon bound states [T2]. They 
are shown to decay algebraically as follows: 

{S^S]) = M^ + Cocos 
ill fl S-+n) = {-iyCi\j\-'/^'^^ + ---, (3) 

n=l n=l 

where M = {Sj) is the uniform magnetization per site. Cm are nonuniversal constants, and k 
is the TL-liquid parameter. On the other hand, correlators of q{< p)-th. multipolar operators 
including transverse spin operator obey an exponential-decay form ~ e"'^'/^ [H] since any 
g-magnon excitation requires a finite excitation energy to break p-magnon bound states. The 
correlation length ^ must be proportional to the inverse of the binding energy A of magnons. 
The absence of the quasi-long-range order in the transverse-spin correlation plays a key role in 
the following discussion of NMR relaxation rate. 

3. NMR relaxation rate 

For spin systems in solids, NMR relaxation time Ti is usually determined through weak hyperfine 
interactions between electron and nuclear spins. Therefore, we can obtain information about 
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electron-spin correlations from measurement of Ti . The formula for Ti |18[ [T9| [20] is given by 



l/Ti 
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where are the hyperfine coupling constants, and S^"{ijj) is the local dynamical structure 
factor 




The strength of depends on kinds of nuclei, crystal structures, and direction of field H. 
The value of u is fixed to a nearly resonating frequency of nuclear-spin precession motion, and it 
is generally much smaller than the energy scale of electron spin systems. Therefore, we take the 
limit u}/{kBT) = uj/3 — > in the following discussions. Utilizing these formulas, we will discuss 
the temperature and field dependence of l/Ti in the octupolar liquid phase. 

4. Temperature and field dependence of 1 /Ti 

From Eqs. ^ and ([5]) and field-theoretical method |18[I19[ [20]. 1/Ti in the multipolar liquid 
phases is obtained as 



Here, /3 = l/(fcBT), u is the excitation velocity, a is lattice spacing, and B{x,y) is the Beta 
function. We note that Eq. ([6]) is valid for the temperature regime /ebT ^ min(n/a. A), where 
the low-energy effective theory of the bound-magnon Bose gas is reliable. If spatial dependence 
of hyperfine couplings is taken account of beyond the formula dH , the prefactors A| in the first 
and second term of Eq. ([6|) are generally modified and become different. 

The form of l/Ti in Eq. ([6]) consists of contributions from the longitudinal spin correlation; 
The first and second terms in Eq. ([6]) originate from the second and third terms in Eq. ([2]), 
respectively. Remarkably, both terms are power-law type functions of temperature. On the 
other hand, since the transverse spin correlation exhibits a gapped behavior, its contribution to 
1/Ti is a thermal-activation form ~ e~^^. Therefore, if we consider the temperature region lower 
than the binding energy A/Ub, the contribution from transverse spin correlation is negligible. 
The disappearance of transverse-spin contribution is a characteristic feature of multipolar TL 
liquids and makes the multipolar liquids distinct from the usual TL liquid in various quasi-lD 
AF magnets as discussed below. 

Parameters k, u, and Cq can be evaluated as a function of field H (or magnetization M) 
and J1/J2 at zero temperature by using density-matrix renormalization group method jl2|, I21|. 
Substituting them into Eq. 1^, we obtain the temperature and field dependence of 1/Ti. In 
particular, the temperature dependence of each term is determined by the TL-liquid parameter 
Hi. For the multipolar phases, k is known to increase almost monotonically with increasing H 
(or M) and approach unity at the saturation [T21II3]- We can thereby find that 1/Ti decreases 
with lowering temperature in the high-field region with k, > 1/2, while it diverges in the low- 
field region with k < 1/2 [14j . This behavior highly contrasts to the usual TL liquid such as 
that of the spin-| AF chain, in which one of the longitudinal- or transverse-spin contributions 
diverges at T — ?• and 1 /Ti is a monotonically increasing function of temperature in sufficiently 
low-temperature regime. Furthermore, for the quadrupolar phase {p = 2), the field dependence 
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Figure 1. Temperature depen- 
dence of NMR relaxation rate 1/Ti 
in the octupolar TL-liquid phase for 
J2/IJ1I = 0.32 (J1/J2 = -3.125). 
We simply set the coupling con- 
stant ^11 = 1. 
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of l/Ti has also been determined rather accurately |15j . It has then been shown that with 
increasing H, 1/Ti decreases for low field regime and increases in the vicinity of the saturation 
field (see Fig. 1 of Ref. [H]). This peculiar non-monotonic field dependence of 1/Ti is expected 
to be a common property of all the multipolar TL-liquid phases in the J1-J2 chain and would 
be also a hallmark of them. 

Let us now examine the temperature and field dependence of 1 /Ti in the octupolar TL-liquid 
phase (p = 3). We evaluate l/Ti combining the formula ([6]) with numerical estimates of k, u, 
and Cq which are obtained in the same manner as Refs. [12^ [21} [22] . The coupling constant Ay is 
simply set to unity. Hereafter, we concentrate on the point of J2/\Ji\ = 0.32 (J1/J2 = —3.125), 
at which the octupolar phase widely exists for 0.07 < M < 1/2 (0.022|Ji| < H < He = 
0.05|Ji|) [12]. We note that in the octupolar phase, the excitation velocity u, i.e., the width of 
low-energy excitation band, is rather small compared to those of quadrupolar phase [15]. For 
instance, we evaluate u/a ~ 0.06|Ji| ~ O.I9J2 for M = 0.2 and u/a ~ 0.01|Ji| ~ O.O3J2 for 
M = 0.4 at the point of J2/\Ji \ = 0.32. This is a natural result because the parameter region of 
the octupolar liquid phase is closer to the ferromagnetic instability point J2/\Ji\ = 0.25 at zero 
field. The magnon binding energy is also evaluated as A ~ 0.025|Ji| ~ 0.078J2 for M = 0.2 
and A ~ 0.041|Ji| ~ O.I3J2 for M = 0.4 at J2/IJ1I = 0.32 [23]. Therefore, when we apply 
the formula ([6]) at this point, temperature must be low enough, say, ksT < 0.01|Ji|. We also 
note that the numerical estimates of Cq in the present calculation are subject to a rather large 
finite-size effect and the evaluated l/Ti is less accurate than that of the quadrupolar liquid in 
Ref. [15]. However, we have checked that the result below is stable and semi-quantitatively 
unchanged against such an error of Cq. 

We show in Fig. [Tj the temperature dependence of 1 /Ti . The result confirms the previous 
prediction that, as temperature decreases, 1/Ti increases divergently for the low-field regime 
(k < 1/2) and decreases for the high-field regime {k > 1/2). Figure [2] shows the field dependence 
of 1/Ti for several fixed temperatures. It is remarkable that, for the sufficiently low temperature 
ksT = 0.001|Ji|, 1/Ti decreases with increasing H (or M) except for the vicinity of the 
saturation field. The contribution from each term in Eq. ^ is plotted in Fig. [3] separately. 
It is clear that the first term in Eq. (jG]) proportional to contributes dominantly to 

a peculiar behavior of 1/Ti, decrease with increasing H, for a wide field range. This feature 
becomes less significant as temperature increases and the field dependence of 1/Ti becomes 
nearly fiat for a rather high temperature, ksT = 0.01|Ji|. The T and H dependence of 1/Ti 
found above is qualitatively the same as that of the quadrupolar TL liquid [15], as expected. Even 
in the neighboring hexadecapolar liquid phase, similar properties must be realized (although the 
formula ([6|) would be valid only in extremely low temperature regime). 

Quasi- ID J1-J2 magnets with high crystallographic symmetry often have principal axes for 
hyperfine coupling tensor. For such a case, we can eliminate all the contributions from the 
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Figure 2. Field dependence of NMR 
relaxation rate 1 /Ti in the octupolar TL- 
liquid phase for J2/\Ji \ = 0.32 (J1/J2 = 
—3.125). Symbols represent the data 
obtained from Eq. ^ with numerical 
values of k, u, and Cq. Lines are guide 
for the eyes. 
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Figure 3. Field dependence of each 
contribution in NMR relaxation rate 
l/Ti [Eq. dlD] in the octupolar TL- 
liquid phase for J2/IJ1I = 0.32 (J1/J2 = 
-3.125) and kBT/\Ji \ = 0.001. 



longitudinal spin correlations, i.e., terms, in 1/Ti, setting the field H parallel to a principal 
axis. As a result, thermal activated contribution ~ e~^^ proportional to A^, which is the 
hallmark of the multipolar liquids, becomes detectable in 1/Ti. In this case, therefore, we 
can easily distinguish the multipolar TL liquid and usual TL liquid phases from temperature 
dependence of 1/Ti. 

5. Discussions 

We have studied the NMR relaxation rate l/Ti in field- induced octupolar TL liquid of the 
spin-i J1-J2 chain. Substituting numerical estimates of parameters u, and Co into the field- 
theoretical result, we have evaluated 1/Ti as a function of temperature and magnetic field. 
We have then found that, similarly in the quadrupolar liquid, 1/Ti in the octupolar TL liquid 
also exhibits characteristic features which are qualitatively different from those in the usual 
TL liquid; In the octupolar liquid, (i) l/Ti decreases (increases) with increasing temperature 
T for the low- (high-)field regime, and (ii) for sufficiently low temperatures, 1/Ti decreases 
with increasing magnetic field H for a wide field range except the vicinity of the saturation 
field. These features can be understood from the absence of the quasi-long-range order in the 
transverse-spin correlation, which is a natural consequence of the formation of magnon bound 
states, and therefore be generally expected in all the multipolar TL liquids in J1-J2 chain. We 
note that the above features of l/Ti in the octupolar phase is observed only in sufficiently low 
temperature regime {ksT < 0.01|Ji| at J2/\Ji\ = 0.32). This is a quantitative difference from 
the case of the quadrupolar phase, which has relatively small | Ji|. 

The result in this paper indicates that NMR relaxation rate 1/Ti could be useful to detect 
a signature of the multipolar TL-liquid phases. However, only the temperature dependence of 
1/Ti cannot be a sufficient evidence for them. For instance, 1/Ti also decreases with lowering 
temperature in high-field spin-gapped states such as magnetization plateaus. Therefore, we 
propose that observations of both characteristic behavior in 1/Ti and gapless behavior of bulk 
quantities (specific heat, susceptibilities, etc.) present an indirect but strong evidence for 
multipolar TL liquids. 

It has been reported [3 El El [HI [9] that three quasi-lD cuprates PbCuS04(OH)2, 



Rb2Cu2Mo30i2, and Li2ZrCu04 have relatively large value of |Ji| (2 < IJ1I/J2 ^ 4). We 
expect that our result will be helpful for searching for the appearance of octupolar phase in 
these compounds under high magnetic field. 
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